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CALIBRATION OF DOWNLINK AND UPLINK CHANNEL 
RESPONSES IN A WIRELESS MIMO COMMUNICATION 

SYSTEM 

BACKGROUND 

Field 

[0001] The present invention relates generally to communication, and more specifically 

to techniques for calibrating downlink and uplink channel responses in a wireless 
multiple-input multiple-output (MIMO) communication system. 

Background 

[0002] A MIMO system employs multiple (Nt) transmit antennas and multiple (Nr) 

receive antennas for data transmission. A MIMO channel formed by the Nt transmit 
antennas and Nr receive antennas may be decomposed into Ns spatial channels, where 
N s < min {N T , N R } . The N s spatial channels may be used to transmit data in parallel 

to achieve higher overall throughput or redundantly to achieve greater reliability. 
[0003] To achieve high performance, it is often necessary to know the response of the 

MIMO channel. For example, an access point may need to know the response of a 
downlink channel in order to perform spatial processing for a downlink transmission to 
a user terminal. In one conventional channel estimation technique, the access point 
transmits a pilot on the downlink, and the user terminal estimates the downlink channel 
response based on the pilot and sends the downlink channel response estimate back to 
the access point. This channel estimation technique utilizes uplink resources and further 
incurs a delay to send back the channel response estimate, both of which are 
undesirable. 

[0004] A TDD system uses a single frequency band for both the downlink and uplink, 

with the downlink being allocated a portion of the time and the uplink being allocated 
the remaining portion of the time. For a TDD system, the downlink and uplink channel 
responses may be assumed to be reciprocal of one another. That is, if H represents a 
channel response matrix from antenna array A to antenna array B, then a reciprocal 
channel implies that the coupling from array B to array A is given by H r , where H r 
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denotes the transpose of H. With a reciprocal channel, the channel response for one 
link (e.g., the downlink) may be estimated based on a pilot received via the other link 
(e.g., the uplink). 

[0005] The access point and user terminal both utilize transmit and receive chains for 

transmission and reception, respectively. A downlink transmission would then observe 
an "effective" downlink channel response that includes the responses of the transmit 
chain at the access point and the receive chain at the user terminal. Correspondingly, an 
uplink transmission would observe an effective uplink channel response that includes 
the responses of the transmit chain at the user terminal and the receive chain at the 
access point. The responses of the transmit and receive chains at the access point are 
typically different from the responses of the transmit and receive chains at the user 
terminal. As a result, the effective downlink and uplink channel responses are typically 
not reciprocal of one another. If the channel response estimate obtained for one link is 
used for spatial processing for the other link, then any difference in the responses of the 
transmit/receive chains at the access point and user terminal would represent error that, 
if not determined and accounted for, may degrade performance. 

[0006] There is, therefore, a need in the art for techniques to calibrate the downlink and 

uplink channel responses in a TDD MIMO system. 



SUMMARY 

[0007] Techniques to calibrate the downlink and uplink channel responses to account 

for differences in the responses of the transmit and receive chains at the access point 
and user terminal are described herein. After calibration, a channel response estimate 
obtained for one link may be used as a channel response estimate for the other link. 
This can simplify channel estimation and spatial processing. 

[0008] The calibration may be separated into two parts - initial calibration and follow- 

on calibration. For the initial calibration, the access point and user terminal transmit 
MIMO pilots (described below) on the downlink and uplink, respectively. The MIMO 
pilots are used to derive "effective" downlink and uplink channel response estimates, 

H dn and H up , which include the responses of the applicable transmit/receive chains. 
The channel estimates and H up are used to derive correction matrices K ap and 
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K ut , which are thereafter used by the access point and user terminal, respectively, to 
account for the responses of their transmit/receive chains, as described below. 
[0009] For the follow-on calibration, one entity (e.g., the access point) transmits a 

MIMO pilot and a steered reference (described below). The other entity (e.g., the user 
terminal) derives (1) an "actual received" transmit matrix V a based on the steered 
reference and (2) a "hypothesized" transmit matrix V hyp based on the MIMO pilot and 
calibration error matrices Q ap and Q ut . The matrices Q ap and Q ut contain guesses or 

estimates of the errors in the correction matrices K ap and K ut , respectively. The 
difference between the transmit matrices V a and V hyp are indicative of the accuracy of 
the estimates of the errors in the correction matrices. The matrices Q ap and Q ut may be 

adjusted based on an adaptive procedure to minimize the error between V a and V hyp . 
Several adaptive procedures to iteratively adjust the matrices Q ap and Q ut are described 

below. The correction matrices K ap and K ut may thereafter be updated by the 

calibration error matrices Q ap and Q ut , respectively. 

[0010] Various aspects and embodiments of the invention are described in further detail 

below. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0011] FIG. 1 shows the transmit and receive portions at the access point and user 

terminal in the TDD MIMO system; 
[0012] FIG. 2 shows the use of the correction matrices at the access point and user 

terminal to account for their transmit/receive chains; 
[0013] FIG. 3 shows a process performed by the access point and user terminal for 

initial calibration, normal operation, and follow-on calibration; 
[0014] FIG. 4 shows a minimum mean square error (MMSE) adaptive procedure; 

[0015] FIG. 5 shows a steepest descent adaptive procedure; and 

[0016] FIG. 6 shows a block diagram of the access point and user terminal. 
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DETAILED DESCRIPTION 

[0017] The word "exemplary" is used herein to mean "serving as an example, instance, 

or illustration." Any embodiment described herein as "exemplary" is not necessarily to 
be construed as preferred or advantageous over other embodiments. 

[0018] The calibration techniques described herein may be used for single-carrier as 

well as multi-carrier TDD MMO systems. For clarity, these techniques are described 
for a single-carrier TDD MIMO system. 

[0019] FIG. 1 shows a block diagram of the transmit and receive portions at an access 

point 110 and a user terminal 150 in a TDD MIMO system 100. For the downlink, at 
access point 110, transmit symbols (denoted by a vector x^) are processed by a 
transmit chain 114 and transmitted from N ap antennas 116 and over a wireless MIMO 
channel having a response of H. At user terminal 150, the N ap downlink signals are 
received by N m antennas 152 and processed by a receive chain 154 to obtain received 
symbols (denoted by a vector ). The processing by transmit chain 1 14 typically 
includes digital-to-analog conversion, amplification, filtering, frequency upconversion, 
and so on. The processing by receive chain 154 typically includes frequency 
downconversion, amplification, filtering, analog-to-digital conversion, and so on. 

[0020] For the uplink, at user terminal 150, transmit symbols (denoted by a vector x up ) 

are processed by a transmit chain 164 and transmitted from N ut antennas 152 and oyer 
the MIMO channel. At access point 110, the N ut uplink signals are received by N ap 
antennas 116 and processed by a receive chain 124 to obtain received symbols (denoted 
by a vector r up ). 

[0021] For the downlink, the receive vector at the user terminal may be expressed as: 

r dn =S ut HT ap x dn , Eq(l) 

where Xdn * s the vector with N ap transmit symbols sent from N ap access point antennas; 

is the vector with N ut received symbols obtained via N ut user terminal antennas; 
T ap is an N ap xN ap diagonal matrix with N ap complex gains for the access point 

transmit chain, one complex gain for each access point antenna; 
R ut is an N ut xN ut diagonal matrix with N ul complex gains for the user terminal 
receive chain, one complex gain for each user terminal antenna; and 
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H is the N ut xN ap channel response matrix for the downlink. 

The responses of the transmit/receive chains and the MIMO channel are typically a 
function of frequency. For simplicity, a flat-fading channel with a flat frequency 
response is assumed. 

[0022] For the uplink, the receive vector at the access point may be expressed as: 

r up =R ap H r T ut x up , Eq(2) 

where x up is the vector with N ut transmit symbols sent from N ut user terminal antennas; 

r up is the vector with N ap received symbols obtained via N ap access point antennas; 
T ut is an N ut x N ut diagonal matrix with N ut complex gains for the user terminal 

transmit chain, one complex gain for each user terminal antenna; 
R ut is an N ap x N ap diagonal matrix with N ap complex gains for the access point 

receive chain, one complex gain for each access point antenna; and 
H r is the N ap x N ut channel response matrix for the uplink. 

[0023] From equations (1) and (2), the "effective" downlink and uplink channel 

responses, and H up , which include the responses of the applicable transmit and 

receive chains, may be expressed as: 

Hdn = R u tHT ap and H up = R ap H r T ut . Eq (3) 

Combining the two equations in equation set (3), the following may be obtained: 

Hip = I ut Rj T ap R ap = Kj H dn K ap or H up = (Kj H dn K ap ) T , Eq (4) 

where K ap = T a ' p R ap and K ul =T^R ut . K ap is an N ap xN ap diagonal matrix for the 
access point and obtained by a ratio of the receive chain response R ap to the transmit 
chain response T ap , where the ratio is taken element-by-element. Similarly, K ut is an 
N ut xN ut diagonal matrix for the user terminal and obtained by a ratio of the receive 
chain response R ut to the transmit chain response T ut . 
[0024] Equation (4) may also be expressed as: 
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H cup = H up K ut = (H dn K ap ) r = H^dn , Eq (5) 

where H cup denotes the calibrated channel response for the uplink; and 
H cdn denotes the calibrated channel response for the downlink. 

The matrices K ap and K ul include values that account for differences in the 

transmit/receive chains at the access point and user terminal. The application of the 
diagonal matrices, K ap and K ut , to the effective downlink and uplink channel 

responses, as shown in equation (5), allows the calibrated channel response for one link 
to be expressed by the calibrated channel response for the other link. 
[0025] Initial calibration may be performed to determine the matrices K ap and K ut . 

Typically, the true channel response H and the transmit/receive chain responses are not 
known nor can they be exactly or easily ascertained. Instead, the effective downlink 
and uplink channel responses, and H up , may be estimated based on MIMO pilots 

sent on the downlink and uplink, respectively. A MIMO pilot is a pilot comprised of N T 
pilot transmissions sent from Nt transmit antennas, where the pilot transmission from 
each transmit antenna is identifiable by the receiving entity. This can be achieved, for 
example, by using a different orthogonal sequence for the pilot transmission from each 
transmit antenna. Estimates of the matrices K ap and K ut (which are called correction 

matrices K ap and K ul ) may then be derived based on the effective downlink and uplink 

channel response estimates, and H up , as described below. The matrices K ap and 

K ut include correction factors that can account for differences in^ the transmit/receive 
chains at the access point and user terminal. 
[0026] FIG- 2 shows the use of the correction matrices K ap and K ut at access point 

110 and user terminal 150. On the downlink, the transmit vector is first multiplied 
with the correction matrix K ap by a unit 112. The subsequent processing by transmit 
chain 114 and receive chain 154 for the downlink is as described in FIG. 1. Similarly, 
on the uplink, the transmit vector x up is first multiplied with the correction matrix K ul 

by a unit 162. The subsequent processing by transmit chain 164 and a receive chain 124 
for the uplink is also as described in FIG. 1. 



* 
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[0027] For a MIMO system, data may be transmitted on Ns eigenmodes of a MIMO 

channel. These eigenmodes may be viewed as orthogonal spatial channels of the 
MIMO channel. The channel response matrix H may be "diagonalized" to obtain the 
Ns eigenmodes of the MIMO channel. This diagonalization may be achieved by 
performing either singular value decomposition of H or eigenvalue decomposition of a 

correlation matrix of H, which is R = H"H, where H" denotes the conjugate 
transpose of H . 

[0028] Table 1 shows the effective and calibrated channel responses for the downlink 

and uplink as well as the singular value decomposition of the calibrated downlink and 
uplink channel response matrices. 



Table 1 - Singular Value Decomposition 





Downlink 


Uplink 


Effective Channel Response 


=R ut HT ap 


Hup = Rap Xut 


Calibrated Channel Response 


Hcdn - HdnKap 


ti _n if 

—cup — up— UI 


Singular Value Decomposition 
of True Channel Response 


Hcdn = Yut? Hap 


H cup =U ap EV u H t 


Singular Value Decomposition 
of Estimated Channel Response 


~ * ~T ~ T 

H cdo =V ut Z U ap 




Unnormalized 
Transmit Matrices 




U ap ? = H CU p ¥-ut 



In Table 1, U ap is an N ap xN ap unitary matrix of left eigenvectors of H C u P > S is an 
N ap xN ut diagonal matrix of singular values of H cup , V ut is an N m xN ut unitary 
matrix of right eigenvectors of H C u P > anc * " * " denotes the complex conjugate. A 
unitary matrix M is characterized by the property M W M = I, where I is the identity 
matrix. The matrices V* t and U* p are also matrices of left and right eigenvectors, 
respectively, of H c <in- The matrices V, V*, V r and V w are different forms of the 
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matrix V. For simplicity, reference to the matrices U ap and V ul in the following 
description may also refer to their other forms. The matrices U ap and V ul (which are 

also called transmit matrices) may be used by the access point and user terminal, 

respectively, for spatial processing and are denoted as such by their subscripts. 
[0029] The singular value decomposition is described in further detail by Gilbert Strang 

entitled "Linear Algebra and Its Applications," Second Edition, Academic Press, 1980, 

which is incorporated herein by reference. 
[0030] In a practical system, the matrices H cdn and H cup are not available. Instead, the 

user terminal can estimate the calibrated downlink channel response based on a M1MO 
pilot transmitted by the access point. The user terminal can then perform singular value 

decomposition of the calibrated downlink channel response estimate, H cdn , to obtain a 

diagonal matrix 2 and a matrix V ut of left eigenvectors of H cdn , where the hat (" A ") 
above each matrix indicates that it is an estimate of the actual matrix. Similarly, the 
access point can estimate the calibrated uplink channel response based on a MIMO pilot 
transmitted by the user terminal. The access point can then perform singular value 

decomposition of the calibrated uplink channel response estimate, H cup , to obtain the 

diagonal matrix £ and a matrix U ap of left eigenvectors of H cup . 

[0031] Because of the reciprocal channel and the calibration, the singular value 

decomposition only needs to be performed by either the user terminal or the access 

point to obtain both matrices V ut and U ap . For clarity, the following description is for 

an implementation whereby the user terminal obtains the calibrated downlink channel 

response estimate, H cdn , performs decomposition of H cdn , uses the matrix V ul for 

spatial processing, and sends the matrix U ap back to the access point using a steered 

reference, as described below. A steered reference (or steered pilot) is a pilot that is 
transmitted from all antennas and on the eigenmodes of the MIMO channel. 
[0032] The user terminal may transmit an uplink steered reference, as follows: 

x U p,m =K ut \ utmPm , Eq(6) 



where p m is a pilot symbol transmitted on eigenmode m for the steered reference; 
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x is a transmit vector for the uplink steered reference for eigenmode m; and 



v ulm is the m-th eigenvector or column of V ut , where V ut = [v uU v ut2 ... v m Nm ] . 



[0033] 



The received uplink steered reference at the access point may be expressed as: 




^HcupiuvnPm + Hup = U ap tY.l l^m P m + & 
= Ma P ,m^mPm + Oup 



Eq(7) 



where 



r upm is a received vector for the uplink steered reference for eigenmode m; 
o m is the m-th diagonal element of % ; and 



u apm is the m-th eigenvector or column of U ap , where U ap = [u ap x u ap2 ... u aPtN-p ] . 

Equation (7) shows that the received uplink steered reference at the access point, in the 
absence of noise, is approximately equal to u apm <7 m /? m - The access point can obtain an 

estimate of the uplink channel response based on the steered reference sent by the user 
terminal using various estimation techniques. 
[0034] In one embodiment, to obtain an estimate of u apm , the received vector r up m is 

first multiplied with the complex conjugate of the pilot symbol, or p* m , and then 
integrated over multiple received steered reference symbols for each eigenmode to 
obtain a vector f upm , which is an estimate of u apm tr m for eigenmode m. Since the 

eigenvectors have unit power, the singular value & m for each eigenmode may be 
estimated based on the received power of the uplink steered reference for that 
eigenmode, which is & m =||r u > m || ). Each of the N ap elements of the estimate of u apm 

is obtained by dividing a corresponding one of the N ap elements of r upm by & m . 
[0035] In another embodiment, an estimate of u apm is obtained based on the received 

vector r up m and using an MMSE technique. Since the pilot symbols p m are known, the 
access point can derive the estimate of u apm such that the mean square error between a 
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received pilot symbol p m (which is obtained after performing the matched filtering on 

the received vector r upm with v utm ) and the transmitted pilot symbol p m is minimized. 

[0036] The access point may perform additional processing on the estimates of u apm , 

for m = 1 .. N s . For example, since these estimates are obtained for one eigenvector at 

a time, the Ns eigenvector estimates may not be orthogonal to one another due to, for 
example, noise in the received steered reference, changes in the MIMO channel 
response, and so on. The access point may then perform Gram-Schmidt 
orthogonalization on the N s eigenvector estimates to obtain orthogonal transmit vectors. 

In any case, the access point obtains a transmit matrix U ap , which is an estimate of U ap , 

which in turn is derived by the user terminal based on H cdn . The access point uses the 

transmit matrix U ap to perform spatial processing for downlink transmission. 

1. Follow-on Calibration 
[0037] The correction matrices K ap and K ut obtained from the initial calibration may 

contain errors due to various sources such as (1) imperfect channel estimates and 
H up used for the initial calibration, (2) changes in the transmit/receive chains at the 
access point and user terminal, and so on. Errors in the correction matrices cause errors 
in both (1) the transmit matrix V ut used by the user terminal for spatial processing and 

derived from H cdn and (2) the transmit matrix U ap used by the access point for spatial 

processing and derived from the uplink steered reference sent using V ut . Improved 
performance may be obtained if the errors in the correction matrices can be estimated 
and removed. 

[0038] The access point and/or the user terminal can perform follow-on calibration to 

estimate the errors in the correction matrices K ap and K ut . For clarity, the following 

description is for follow-on calibration by the user terminal. For follow-on calibration 
by the user terminal, the access point transmits a MIMO pilot on the downlink using the 

correction matrix K ap and also transmits a steered reference on the downlink using the 
transmit matrix U nn and the correction matrix K nn . The downlink steered reference 
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may be expressed as: = K ap u apm /7 m , where t) ap =[fi ap l u ap2 ... u ap , N J. The user 

terminal can obtain an estimate of V ut Z r based on the received downlink steered 
reference, similar to that described above for the uplink steered reference. For 
simplicity, the estimate of V ut S r derived from the downlink steered reference is called 
an "actual received" transmit matrix V a , which is an unnormalized matrix that includes 
an estimate of V ut as well as an estimate of 2 . (For clarity, the "~" above a matrix 
indicates that it is an unnormalized matrix.) The user terminal also obtains another 
version of H cdn from the MIMO pilot sent by the access point. 
[0039] The errors in the correction matrices K ap and K ut may be represented by 

diagonal calibration error matrices and , respectively. The correction matrices 

K ap and K ut may then be expressed as: 

K ap = K ap Q' ap and K ut = K ut Q' m . Eq (8) 

If the errors in the correction matrices are small, then the diagonal elements of Q^ ap and 
Q' m are complex values close to 1 + jO . The calibrated downlink channel response 
estimate, H cdn , may then be expressed as: 

H cdn = H^K.p = IL K ap Q' ap , Eq (9) 

or H^Q J^H^K.p . 

[0040] The matrices Q' and contain the "true" errors in K ap and K ul , 

respectively. A guess or estimate of Q[ ap and may be denoted as Q ap and Q m , 
respectively. A "hypothesized" downlink channel may be defined as: 

H hyp =H cdn Q a -; . Eq(10) 
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A hypothesized downlink channel is a guess of H^K^ and is derived under the 

assumption that the error in the applied correct correction matrix K op is Q ap . If Q ap is 

a perfect guess of Q[ ap in equation (10) and is a perfect estimate of in 

equation (9), then H hyp = H cdn and = H cup . 
[0041] The spatial processing at the access point may be expressed as: 



£ ap | = H cup V ut = H up K ut V ul = H up K ut Q' u , V ul , Eq (1 1) 

where V ut is obtained from the singular value decomposition of H cdn , which is 
obtained from the downlink MIMO pilot. The user terminal does not have the value for 
Q' ap , but only its guess Q ap . The user terminal thus computes an unnormalized 

transmit matrix that hypothetically would have been obtained by the access point if 
the calibration error matrices were Q and Q , as follows: 

— ap _ut ' 

£ ra =H: yp Q ui V ut =(H e<ta Q;;) r Q ut V ut . Eq (12) 

Equation (12) is equal to equation (11) if Q ap is a perfect guess of Q[ ap and Q ut is a 

perfect guess of Q' . If this is the case, then H hyp = H^K^ . 

[0042] The user terminal then performs processing on in the same manner that the 

access point would have performed on a received uplink steered reference and obtains a 
"generated" transmit matrix U g , which is a normalized transmit matrix that resembles 

U ap . For example, the access point may perform Gram-Schmidt orthogonalization of 

the received eigenvectors u apm in order to improve the performance of its transmit 
steering vectors. In this case, the user terminal would perform the same 
orthogonalization on the eigenvectors in . The user terminal simply emulates the 
processing that is normally performed by both the access point and the user terminal, 
albeit under an assumption of calibration errors represented by Q ap and Q m . The 
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matrix U g would have been used by the access point to transmit the downlink steered 

reference and for spatial processing of downlink transmission. 
[0043] The spatial processing at the user terminal may be expressed as: 

Y,g = HcanUg = H dn K ap U g = H dn K ap Q' ap U g . Eq (13) 

Again, the user terminal does not have Q' ap , but only its guess Q ap . The user terminal 
thus computes a hypothesized transmit matrix V hyp for itself as follows: 

V hyp = H hyp Q ap U g . Eq(14) 

Equation (14) is equal to equation (13) if H hyp is a perfect guess of H^K^ and Q ap is 
a perfect guess of Q^ ap . The matrix V hyp is an unnormalized matrix that includes a user 
terminal transmit matrix V g (which corresponds to the access point transmit matrix 
U g ) as well as a diagonal matrix Z g (which resembles 2). The matrix V hyp is 
hypothesized to have been received by the user terminal with (1) the user terminal 
transmitting an uplink steered reference using V ut , (2) the access point performing its 
normal processing on the received uplink steered reference to derive its transmit matrix 
U g , (3) the access point transmitting a downlink steered reference using U g , (4) the 

correction matrices K ap and K ut having the errors indicated by the matrices Q ap and 

Q m , respectively, and (5) assuming no channel estimation error in H cdn from the 
downlink MIMO pilot. 

[0044] Equations (12) and (14) are correct if the calibration error matrices Q ap and Q m 

correctly indicate the true errors in the correction matrices K ap and K ul , respectively. 
The difference between the actual received transmit matrix V a obtained from the 
downlink steered reference and the hypothesized transmit matrix V hyp obtained from 
the downlink MIMO pilot may be computed as follows: 

E = V a -V hyp) Eq(15) 
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where E is an N ut xN ap matrix of errors between V a and Y hyp . The error matrix E 

gives an indication of the accuracy of the guess for Q ap and Q ut . Various adaptive 

procedures may be used to adjust the matrices Q ap and Q m to drive the error matrix E 

toward zero. Such adaptive procedures include an MMSE adaptive procedure and a 
steepest descent adaptive procedure. The diagonal elements of Q ap and Q m may be 

initialized to 1 + jO for the adaptive procedures. 
[0045] For the MMSE adaptive procedure, approximate partial derivatives of the 

elements of E are computed with respect to the elements of Q ap and Q m . If the "lead" 

element (which is the upper leftmost element) of K ap is set to 1+ jO by the initial 

calibration, then this element does not need to be adjusted. In addition, the error matrix 
E is not affected by the magnitude of the elements in Q ut . Thus, Q ut may be 

normalized, for example, by defining the real component of the lead element of Q ut to 

be 1.0. Furthermore, an eigenvector may be scaled by any unit-magnitude complex 
number (or rotated by any phase) without affecting its performance. Thus, a set of 
phases may be selected to make V hyp close to V a without any loss in generality. This 
property allows for scaling of Q m by an arbitrary unit-magnitude factor, so the 
imaginary component of the lead element of Q m may be defined to be 0.0. 
[0046] The MMSE adaptive procedure may be performed as follows. Let q be a real 

vector of length 2(N ap + N ut - 2) and made up of the real and imaginary components of 
the diagonal elements of Q ap and Q ut , except for the lead elements which are set to 1.0. 
The vector q may be defined as: 

?2o-2 )+ i = Re{Q ap (iJ)} , for * = 2 ... N ap , 

?2(,-i) =Jm iQop (M)} , for i = 2 ... N ap , 

?2(/ + n v -2h =Re{Q tt ,(M)} , for i = 2 ... N ut , and 

420+^-2) =Im{Q M/ (/,i)} , for i = 2 ... N ut , 



where q i is the z-th element of q ; 
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Q ap (i,i) is the i-th diagonal element of Q ap ; and 
Q ut (iJ) is the i-th diagonal element of Q m . 

The odd-indexed elements of q are for the real component of the diagonal elements of 
Q ap and Q ut , and the even-indexed elements of q are for the imaginary component of 
the diagonal elements of Q ap and Q m . The first 2N ap - 2 elements of q are for the 
N ap -1 diagonal elements other than the lead element of Q ap , and the last 2N ut -2 
elements of q are for the N ut - 1 diagonal elements other than the lead element of Q m . 
[0047] Let e be a real vector of length 2N ap N ut and made up of the real and imaginary 

components of the elements of E . The vector e may be defined as: 

^()-i)N tt )-i =Re { £ ftj)} > for i = 1 ... N ut and j = l ... N ap , and 
e 2UHj-i)*«) =kn{E(iJ)} , for i = l ... N ut and j = l ... N ap , 

where e i is the *-th element of e ; and 

£(/, j) is the element in the /-th row and j-th column of E . 

The odd-indexed elements of e are for the real component of the elements of E , and 
the even-indexed elements of e are for the imaginary component of the elements of E . 
The error vector e can be obtained by evaluating equations (10), (12), (14), and (15) 
with the vector q . 

[0048] For the MMSE adaptive procedure, the partial derivative of a real or imaginary 

component of an element in E with respect to a real or imaginary component of an 
element in Q ap or Q m may be generated by perturbing the component of the element in 

Q ap or Q m and computing the function defined by equations (10), (12), (14), and (15). 

As part of the computation for Y hyp , a single e jx term may be selected such that 

I Y a ~e JX ' Yhyp P is minimized. This is done to normalize the phase of the lead element 
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[0049] Approximate partial derivatives of the elements of e with respect to the 

elements of q may be expressed as: 

a ^ j gj .(q + A,)-gy(q) for «' = 1 - 2(N ap + N M -2) 

A. . = — ^ = - J —= , Eq (16) 

" dq t S and y = l ... 2N ap N ut 

where A ; is a vector of length 2(N ap + N ut -2) and containing a small real value of 5 

for the j-th element and zeros elsewhere; and 
Ajj is the approximate partial derivative of the 7-th element of e with respect to 

the i-th element of q . 

The approximate partial derivative A, . may be obtained as follows. A vector q. is first 
computed as q = q + A f . The function defined by equations (10), (12), and (14) is then 
evaluated for q. (which contains Q and Q uM ) to obtain a new (or "revised") 
hypothesized transmit matrix V hyp . Y hyp f is then subtracted from V a to obtain a new 
error matrix E r = V a - V hyp / , which is used to form a new error vector e, . The 7-th 
element of e , which is denoted as <?;(q) in equation (16), is then subtracted from the 7- 
th element of e, , which is denoted as e .(q + A f .) in equation (16). The result of the 
subtraction is divided by Sto obtain A j t . . 

[0050] The computation for equations (10), (12), (14), and (15) is performed for each of 

the 2(N ap + N ut - 2) elements of q to obtain a corresponding new error vector e, . For 

each new error vector e, , the 2N ap N ut elements of e are subtracted from the 2N ap N m 

elements of e r , on an element-by-element basis, to obtain 2N ap N ut approximate partial 

derivatives of the 2N ap N ut elements of e with respect to the /-th element of q. A 

matrix A of dimension 2N ap N ut by 2(N ap + N ul - 2) may be formed with all of the 

partial derivatives for all of the elements of e and q. Each column of A contains 

2N ap N ut approximate partial derivatives for the 2N ap N ul elements of e with respect to 

one element of q. The 2(N ap + N ut -2) columns of A are for the 2(N ap +N ut -2) 

elements q. 
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[0051] If the relationships in equations (10), (12), (14), and (15) are approximately 

linear, then an estimate of the difference between the guess of the calibration errors in q 

and the actual calibration errors may be expressed as: 

y=A" 1 e , Eq(17) 

where y is an update vector for the estimated difference between q and the actual 
calibration errors. The update vector y has the same format and dimension as the 
vector q , which is a real vector made up of the real and imaginary components of the 
diagonal elements of Q ap and Q ut other than the lead elements. 

[0052] If A is not a square matrix, which is typically the case, then the simple matrix 

inverse does not exist. The Moore-Penrose pseudo-inverse of A may then be used for 
equation (17). This pseudo-inverse is simply a matrix that satisfies the equations 
AA" 1 A = A and A" 1 AA -1 = A" 1 . The pseudo-inverse may be generated by performing 
singular value decomposition of A , which is A = U„ DV " , and computing the pseudo- 
inverse as A -1 = U"D _1 V a , where D" 1 is a diagonal matrix made up of the inverses of 
the corresponding non-zero diagonal elements of D . 

[0053] The matrix A of partial derivatives is computed under an assumption that the 

function defined by equations (10) through (13) is approximately linear for calibration 
errors of the size being evaluated. Since the linearity assumption is not completely 
accurate, the procedure may be iterated multiple times to determine the correct 
calibration errors. For some cases, the procedure does not converge. However, 
convergence can generally be achieved by simply selecting a different initial guess for 
the calibration errors. If convergence is not obtained, the user terminal can also obtain 

another version of V a and H cdn based on another estimate of the downlink steered 
reference and downlink MIMO pilot and perform the MMSE adaptive procedure using 
these new matrices. 

[0054] If equations (10), (12), (14), and (15) are linear, then y +q would minimize the 

mean square of the elements of e. However, since these equations are not linear, q is 
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replaced with y + q and the procedure is repeated. The calibration error vector may 
then be updated as follows: 

— mmse — — mmse 

where n is an index for the iteration number; 

q (n) is the calibration error vector for the n-th iteration; 

—mmse * 

y(n) is the update vector obtained for the n-th iteration; and 

q (n + 1) is the calibration error vector for the (n + 1) -th iteration. 

— mmse 

[0055] The computation described above may be repeated for any number of iterations. 

Each iteration uses the updated calibration error vector q (n + 1) obtained from the 

* -±mmse 

prior iteration. The procedure can terminate when the update vector y(n) is sufficiently 

small. For example, the termination condition may be || y(n) || 2 = y < y thl , where y is the 

sum of the squares of the magnitude of the elements of y(n) and y thl is a threshold 

value. As another example, the termination condition may be y i < y th2 , for all i, where 

y ( is the i-th element of y(n) and y th2 is another threshold value. After all iterations 

have been completed, the matrices for the final estimates of calibration errors are 
denoted as Q fi , and Q a , . 

_ap,y?nn/ —ut, final 

[0056] For the steepest descent adaptive procedure, an aggregate error may be defined as: 

z=||E|| 2 = ||V a -V hyp || 2 . Eq(19) 

The aggregate error, z, is obtained by summing the squares of the magnitude of the 
elements of E . The partial derivatives of z with respect to the elements of q may be 

computed as follows: 

dz z(Q + A,)-z(q) 
g| =jp-5 a =j ~ , for i = l... 2(N ap + N ut -2), Eq (20) 

where g { is the approximate partial derivative of z with respect to the i-th element of q 
and A, is a vector of length 2(N ap +N ut -2) and containing a small real value of <?for 
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the i-th element and zeros elsewhere. The approximate partial derivative g. may be 
obtained as follows. A vector q is first computed as q = q + A, . The function defined 
by equations (10), (12), (14), and (15) is evaluated for q to obtain a new error vector 
e, . The aggregate error z i is then computed for the new error vector e, , as shown in 
equation (19). The aggregate error z obtained with q, which is denoted as z(q) in 
equation (20), is then subtracted from the aggregate error z f obtained with q , which is 
denoted as z(q + A f ) in equation (20). The result of the subtraction is divided by 8\o 
obtain g i . The computation is repeated for each of the 2(N ap + N ut - 2) elements of q . 
A vector g of dimension 2(N ap + N ut - 2) may be formed with the approximate partial 
derivates obtained for the 2(N ap + N ut -2) elements of q. Each element of g is the 
slope of the aggregate error evaluated at a corresponding element of q . 
[0057] The calibration error vector may then be updated as follows: 

where g(n) is the slope vector obtained for the n-th iteration, and Q jrf (rc) and Q sd (n + 1) 

are the calibration error vectors for the rc-th and (n + 1) -th iterations, respectively, for 
the steepest descent procedure. The computation described above may be repeated for 
any number of iterations. Each iteration uses the updated calibration error vector 
q ^(nH-l) obtained from the prior iteration. The procedure can terminate when the 

aggregate error z is sufficiently small, e.g., less then a z th threshold value. 

[0058] Two adaptive procedures have been described above for deriving estimates of 

the actual calibration errors. Other adaptive and non-adaptive procedures may also be 
used, and this is within the scope of the invention. 

[0059] The user terminal can update its correction matrix to account for the calibration 

errors, as follows: 

K ut ,_=K ut Q^ nfl ,. Eq(22) 
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The user terminal uses the new correction matrix K ut new , instead of the prior correction 
matrix K ut , for spatial processing for uplink transmission, as shown in FIG. 2. The user 
terminal may send the calibration error matrix Q apy?nn/ to the access point, which may 

then update its correction matrix as K an „„ w = K nn Q"' ,. To reduce the amount of 
signaling, the user terminal may only send back the calibration error matrix Q npfintd if 

the matrix meets some predetermined threshold. 
[0060] FIG. 3 shows a process 300 performed by the access point and user terminal for 

initial calibration, normal operation, and follow-on calibration. The access point and 
user terminal perform initial calibration to calibrate their transmit and receive chains 

and derive correction matrices K ap and K ut (block 310). The initial calibration is 

described below. 

[0061] Thereafter, for normal operation, the access point transmits a downlink MIMO 

pilot using its correction matrix K ap (block 322). The user terminal obtains a calibrated 

downlink channel response estimate, H cdn , based on the downlink MIMO pilot (block 
324) and performs singular value decomposition of H cdn to obtain its transmit matrix 
V ut (block 326). The user terminal then transmits an uplink steered reference using 
V ut and K ut , as shown in equation (6) (block 328). The access point receives the 

uplink steered reference and derives its transmit matrix U ap , as described above (block 

330). The access point and user terminal use the transmit matrices U ap and V ut , 

respectively, for spatial processing. 
[0062] For the follow-on calibration, the access point transmits a downlink steered 

reference using U ap and K ap , and further transmits a downlink MIMO pilot using K ap 

(block 342). The user terminal derives the actual unnormalized transmit matrix V a 
based on the downlink steered reference, as described above (block 344). The user 
terminal also computes the unnormalized transmit matrix based on its transmit 
matrix V ut , the calibrated downlink channel response estimate H cdn , and the calibration 
error matrices Q a and Q m , as shown in equations (10) and (12), or 
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£rx =(HcdnQ~p) 7 Q ut Yui (block 346). The user terminal further processes in the 

same manner as would have been performed by the access point (e.g., perform 
orthogonalization) to obtain the normalized transmit matrix U g (block 348). The user 

terminal then computes the hypothesized unnormalized transmit matrix V hyp based on 

the transmit matrix U g and the calibrated downlink channel response estimate H cdn , as 

shown in equations (10) and (14), which is V hyp = H cdn U g (block 350). The matrix 

V hyp is the unnormalized transmit matrix that the user terminal would have received if 

the access point transmits a downlink steered reference using U g . The user terminal 

then revises the calibration error matrices Q and Q based on the transmit matrices 

V a and V hyp (block 352). Blocks 346 through 352 may be performed using an 

adaptive procedure. The user terminal may thereafter update its correction matrix K ut 
with the calibration error matrix Q m (block 354), and the access point may thereafter 

update its correction matrix K ap with the calibration error matrix Q ap (block 356). 

[0063] FIG. 4 shows an MMSE adaptive procedure 400, which may be used for blocks 

346 through 352 in FIG. 3. The hypothesized transmit matrix V hyp is first computed 

based on H cdn and Q ap and Q ut (block 410). Block 410 corresponds to blocks 346 

through 350 in FIG. 3. The error matrix E is next computed as the difference between 
the transmit matrices V a and V hyp , as shown in equation (15) (block 412). Partial 

derivatives for each of the elements in the error matrix E with respect to each of 
selected ones of the elements (e.g., all diagonal elements except for the lead elements) 
in the calibration error matrices Q ap and Q m are then derived, as shown in equation 

(16) (block 414). The matrix E and the matrices Q ap and Q m may be placed in the 

form of vectors for ease of computation, as described above. The partial derivatives 
may be derived separately for the real and imaginary components of the elements in the 
matrices, as also described above. The update vector y is then computed based on the 

matrix A of partial derivatives and the error matrix E, as shown in equation (17) 
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(block 416). The calibration error matrices Q ap and Q m are then updated with the 

update vector y , as shown in equation (18) (block 418). A determination is next made 

whether or not the update vector y satisfies a termination condition (block 420). If the 

answer is 'yes', then process 400 terminates. Otherwise, the process returns to block 
410 and performs another iteration. 
[0064] FIG. 5 shows a steepest descent adaptive procedure 500, which may also be 

used for blocks 346 through 352 in FIG. 3. The hypothesized transmit matrix V hyp is 

first computed based on H cdn and Q ap and Q ut (block 510). The aggregate error z is 

next computed as z - ||V a -V hyp || 2 , as shown in equation (19) (block 512). Partial 

derivatives for the aggregate error with respect to each of selected ones of the elements 
in the calibration error matrices Q ap and Q m are then derived, as shown in equation 

(20) (block 514). The matrices Q ap and Q m may be placed in the form of a vector, and 

the partial derivatives may be derived separately for the real and imaginary components 
of the elements in the matrices. The calibration error matrices Q and Q are then 

— ap — m 

updated with the partial derivatives, as shown in equation (21) (block 516). A 
determination is next made whether or not the aggregate error z satisfies a termination 
condition (block 518). If the answer is 'yes', then process 500 terminates. Otherwise, 
the process returns to block 510 and performs another iteration. 
[0065] In the above description, the user terminal estimates the calibration errors in 

both the correction matrices K ap and K ut . To simplify the follow-on calibration, the 

user terminal can assume that the correction matrix K ap contains no errors and only 

estimate the errors in the correction matrix K ut . This is equivalent to setting the 
calibration error matrix Q ap to the identity matrix. By omitting Q ap , the dimensions of 

the vectors q , y , and g and the matrix A are reduced, which may then greatly reduce 
the computation. 

[0066] For clarity, the description above is for the case in which the user terminal 

performs follow-on calibration. The access point may also perform the follow-on 
calibration. In this case, the access point and user terminal switches role in FIG. 3. The 



040159 

EL977102632US 

23 

user terminal would then transmit an uplink steered reference and an uplink MIMO 
pilot, and the access point would perform the computation to derive Q ap and Q m . 

[0067] Also for clarity, the follow-on calibration is described for a single-carrier MIMO 

system. The follow-on calibration may also be performed for a multi-carrier MIMO 
system, which may utilize orthogonal frequency division multiplexing (OFDM) or some 
other multi-carrier modulation technique. OFDM effectively partitions the overall 
system bandwidth into multiple (N F ) orthogonal subbands, which are also called tones, 
subcarriers, bins, and frequency channels. With OFDM, each subband is associated 
with a respective subcarrier that may be modulated with data. For a MIMO system that 
utilizes OFDM (an MIMO-OFDM system), the computation described above may be 
performed for each of multiple subbands. Since a high degree of correlation may exist 
between nearby subbands, the calibration may be performed in a manner to take 
advantage of this correlation, for example, to reduce the number of subbands to perform 
follow-on calibration, to speed up convergence, and so on. 

2. Initial Calibration 

[0068] For the initial calibration to derive the correction matrices K ap and K ut , one 

entity (either the user terminal or access point) obtains both the effective downlink 
channel response estimate, H^, and the effective uplink channel response estimate, 
H up . The channel estimates and H up may be obtained based on downlink and 
uplink MIMO pilots, respectively. The correction matrices may be computed from 
and H up and using matrix-ratio computation or MMSE computation. 

[0069] For the matrix-ratio computation, an N ul x N ap matrix C is first computed as: 

C = ^ , Eq(23) 

where the ratio is taken element-by-element. 
[0070] The diagonal elements in the correction matrix K ap for the access point are set 

equal to the mean of the normalized rows of C . Each row of C is first normalized by 
scaling each of the N ap elements in that row with the first element in the row. The mean 
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of the normalized rows (denoted by a vector c row ) is then computed as the sum of the 
N m normalized rows divided by N ut . The N ap diagonal elements of K ap are then set 
equal to the N ap elements of c row . Because of the normalization, the lead element of 
K ap is equal to unity. 

[0071] The diagonal elements in the correction matrix K ut for the user terminal are set 

equal to the mean of the inverses of the normalized columns of C . The j-th column of 
C , for j* = l... N ap , is first normalized by scaling each element in that column with the 

y-th diagonal element of K ap . The mean of the inverses of the normalized columns 

(denoted by a vector c col ) is then computed by (1) taking the inverse of each normalized 
column, where the inversion is performed element-wise, (2) summing the N ap inverse 
normalized columns, and (3) dividing each element in the resultant column by N ap to 

obtain c col . The N ut diagonal elements of K ut are set equal to the N ut elements of c co , . 
[0072] For the MMSE computation, the correction matrices K ap and K ut are derived 

from the effective downlink and uplink channel response estimates, and H up , such 

that the mean square error (MSE) between the calibrated downlink and uplink channel 
responses is minimized. This condition may be expressed as: 

mm ICH^K^) 7 -H up K ut | 2 , Eq (24) 

which may also be written as: 



min 



KapHdn ~ Hup Kut 



where K ap = K ap since K ap is a diagonal matrix. 

[0073] Equation (24) is subject to the constraint that the lead element of K ap is set 

equal to unity. Without this constraint, the trivial solution would be obtained with all 
elements of the matrices K ap and K ut set equal to zero. In equation (24), a matrix Y is 

first obtained as Y = K ap H^ -H up K ut . The square of the absolute value is next 



040159 

EL977102632US 

25 

obtained for each of the N ap N ut elements of Y . The mean square error (or the square 
error since a divide by N ap N ut is omitted) is then equal to the sum of all squared values. 
[0074] The MMSE computation is performed as follows. For simplicity, the elements 

of are denoted as {a jy }, the elements of H up are denoted as {b^} , the diagonal 

elements of K ap are denoted as {w f } , and the diagonal elements of K ul are denoted as 
{Vj} 9 where i = l ... N ap and j- \ ... N ut . The mean square error may be rewritten 
from equation (24), as follows: 

MSE = ]T S a v»i ~ V; I 2 ■ E q (26) 



again subject to the constraint u x = 1 . The minimum mean square error may be 
obtained by taking the partial derivatives of equation (25) with respect to u and v and 
setting the partial derivatives to zero. The results of these operations are the following 
equations: 

&w - V>) ■ a l = 0 ' for i = 2 - N a P > and E q ( 26a > 

^(fl^-Vy)-*^ 0 - for 7 = 1 ... N m . Eq(26b) 

In equation (26a), m, = 1 so there is no partial derivative for this case, and the index i 
runs from 2 through N ap . 
[0075] The set of N ap + N ut - 1 equations in equation sets (26a) and (26b) may be more 

conveniently expressed in matrix form, as follows: 

Bk=? , Eq(27) 

where 
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0 



Zh 
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0 -*21«21 



1*1 



0 



0 



1=1 



0 



Nap 

i=i 



k = 



" u 2 ' 




0 


u 3 




0 


v l 


and z = 


0 

a n b n 


V 2 




a i2 b n 






_ a iN u ,*lN ul _ 



The matrix B includes N ap +N ut -1 rows, with the first N — 1 rows 
corresponding to the N ap - 1 equations from equation set (26a) and the last N ut rows 
corresponding to the N ut equations from equation set (26b). The elements of the matrix 
B and the vector z may be obtained based on the elements of and H up . The 

diagonal elements of K ap and K ul are included in the vector k , which may be obtained 

as: 



k = B ! z . 



Eq (28) 
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The results of the MMSE computation are correction matrices K ap and K ul that 

minimize the mean square error in the calibrated downlink and uplink channel 
responses, as shown in equation (24). 

3. Spatial Processing 

[00771 Table 2 summarizes the spatial processing performed by the user terminal and 

access point for data transmission and reception on the eigenmodes of the MIMO 
channel. 

Table 2 





Uplink 


Downlink 


User Terminal 


Transmit : 
x =K V s 

—up — ut — ut -up 


Receive : 

A-l AT 


Access Point 


Receive : 
i U p=^ U ap r up 


Transmit : 

2£dn = K.apHap§dn 



For the embodiment shown in FIG. 2 and Table 2, the correction matrices K ap and K ut 
are applied on the transmit side at the access point and user terminal, respectively. 

4. System 

[0078] FIG. 6 shows a block diagram of an embodiment of access point 110 and user 

terminal 150 in TDD MIMO system 100. On the downlink, at access point 110, a 
transmit (TX) data processor 610 receives traffic data (i.e., information bits) from a data 
source 608 and signaling and other data from a controller 630. TX data processor 610 
formats, codes, interleaves, and modulates (or symbol maps) the different types of data 
and provides data symbols. As used herein, a "data symbol" is a modulation symbol for 
data, and a "pilot symbol" is a modulation symbol for pilot. The pilot symbols are 
known a priori by both the access point and user terminal. A TX spatial processor 620 
receives the data symbols from TX data processor 610, performs spatial processing on 
the data symbols, multiplexes in pilot symbols as appropriate (e.g., for channel 
estimation, calibration, and so on), and provides N ap streams of transmit symbols to N ap 
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modulators (MOD) 622a through 622ap. Each modulator 622 receives and processes a 
respective transmit symbol stream to obtain a corresponding stream of OFDM symbols, 
which is further processed by a transmit chain within the modulator to obtain a 
corresponding downlink modulated signal. N ap downlink modulated signals from 
modulator 622a through 622ap are then transmitted from N ap antennas 624a through 
624ap, respectively. 

[0079] At user terminal 150, N ut antennas 652a through 652ut receive the transmitted 

downlink modulated signals, and each antenna provides a received signal to a respective 
demodulator (DEMOD) 654. Each demodulator 654 (which includes a receive chain) 
performs processing complementary to that performed at modulator 622 and provides 
received symbols. A receive (RX) spatial processor 660 then performs receiver spatial 
processing on the received symbols from all N ut demodulators 654 to obtain detected 
symbols, which are estimates of the data symbols sent by the access point. For follow- 
on calibration, RX spatial processor 660 provides (1) a calibrated downlink channel 

response estimate, H cdn , obtained from a downlink MIMO pilot transmitted by the 
access point and (2) received symbols for a downlink steered reference transmitted by 
the access point. An RX data processor 670 processes (e.g., symbol demaps, 
deinterleaves, and decodes) the detected symbols and provides decoded data. The 
decoded data may include recovered traffic data, signaling, and so on, which are 
provided to a data sink 672 for storage and/or a controller 680 for further processing. 
[0080] The processing for the uplink may be the same or different from the processing 

for the downlink. Data and signaling are processed (e.g., coded, interleaved, and 
modulated) by a TX data processor 690, and further spatially processed and multiplexed 
with pilot symbols by TX spatial processor 692 to obtain transmit symbols. The 
transmit symbols are further processed by modulators 654a through 654ut to obtain N ut 
uplink modulated signals, which are then transmitted via N ut antennas 652a through 

652ut to the access point. User terminal 150 sends back the correction K ap for the 
initial calibration and may send back the calibration error matrix Q , , for the follow- 

J — ap, final 

on calibration, for the implementation described above. At access point 110, the uplink 
modulated signals are received by antennas 624, demodulated by demodulators 622, and 
processed by an RX spatial processor 640 and an RX data processor 642 in a 
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complementary to that performed at the user terminal. RX data processor 642 provides 
the matrices K an and Q , to controller 630. 

[0081] For the initial calibration, the access point and user terminal transmit MIMO 

pilots on the downlink and uplink, respectively. Each entity derives the effective 
channel response estimate for its link. One entity (e.g., the access point) sends the 
channel estimate to the other entity (e.g., the user terminal) for computation of the 

correction matrices K ap and K ut for both entities. The entity that derives the correction 

matrices uses its correction matrix and sends the other correction matrix back to the 
other entity. For the follow-on calibration, one entity (e.g., the access point) transmits 
both the steered reference and MIMO pilot. The other derives the calibration error 
matrices Q apy?mj/ and Q utfinal for both entities based on the received pilots, as described 

above. The entity that derives the calibration error matrices uses its calibration error 
matrix and may send the other calibration error matrix back to the other entity (e.g., if 
the errors are sufficiently large). 
[0082] Controllers 630 and 680 control the operation of various processing units at the 

access point and user terminal, respectively. Controllers 630 and/or 680 may also 
perform processing for the initial and follow-on calibration (e.g., the computation for 

the correction matrices K aD and K m and the calibration error matrices Q a , 
andQ utfinal ). Memory units 632 and 682 store data and program codes used by 

controllers 630 and 680, respectively. 

[0083] The calibration techniques described herein may be implemented by various 

means. For example, these techniques may be implemented in hardware, software, or a 
combination thereof. For a hardware implementation, the processing units used to 
perform the initial and/or follow-on calibration may be implemented within one or more 
application specific integrated circuits (ASICs), digital signal processors (DSPs), digital 
signal processing devices (DSPDs), programmable logic devices (PLDs), field 
programmable gate arrays (FPGAs), processors, controllers, micro-controllers, 
microprocessors, other electronic units designed to perform the functions described 
herein, or a combination thereof. 

[0084] For a software implementation, the calibration techniques may be implemented 

with modules (e.g., procedures, functions, and so on) that perform the functions 
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described herein. The software codes may be stored in a memory unit (e.g., memory 
unit 632 or 682 in FIG. 6) and executed by a processor (e.g., controller 630 or 680). 
The memory unit may be implemented within the processor or external to the processor, 
in which case it can be communicatively coupled to the processor via various means as 
is known in the art. 

[0085] Headings are included herein for reference and to aid in locating certain 

sections. These headings are not intended to limit the scope of the concepts described 
therein under, and these concepts may have applicability in other sections throughout 
the entire specification. 

[0086] The previous description of the disclosed embodiments is provided to enable any 

person skilled in the art to make or use the present invention. Various modifications to 
these embodiments will be readily apparent to those skilled in the art, and the generic 
principles defined herein may be applied to other embodiments without departing from 
the spirit or scope of the invention. Thus, the present invention is not intended to be 
limited to the embodiments shown herein but is to be accorded the widest scope 
consistent with the principles and novel features disclosed herein. 

[0087] WHAT IS CLAIMED IS: 



